Patterns of item nonresponse behaviour to survey questionnaires are systematic and associated with genetic loci.

Response to survey questionnaires is vital for social and behavioural research, and most analyses assume full and accurate response by participants. However, nonresponse is common and impedes proper interpretation and generalizability of results. We examined item nonresponse behaviour across 109 questionnaire items in the UK Biobank (N = 360,628). Phenotypic factor scores for two participant-selected nonresponse answers, 'Prefer not to answer' (PNA) and 'I don't know' (IDK), each predicted participant nonresponse in follow-up surveys (incremental pseudo-R2 = 0.056), even when controlling for education and self-reported health (incremental pseudo-R2 = 0.046). After performing genome-wide association studies of our factors, PNA and IDK were highly genetically correlated with one another (rg = 0.73 (s.e. = 0.03)) and with education (rg,PNA = -0.51 (s.e. = 0.03); rg,IDK = -0.38 (s.e. = 0.02)), health (rg,PNA = 0.51 (s.e. = 0.03); rg,IDK = 0.49 (s.e. = 0.02)) and income (rg,PNA = -0.57 (s.e. = 0.04); rg,IDK = -0.46 (s.e. = 0.02)), with additional unique genetic associations observed for both PNA and IDK (P < 5 × 10-8). We discuss how these associations may bias studies of traits correlated with item nonresponse and demonstrate how this bias may substantially affect genome-wide association studies. While the UK Biobank data are deidentified, we further protected participant privacy by avoiding exploring non-response behaviour to single questions, assuring that no information can be used to associate results with any particular respondents.

Cross-disorder genetic analysis of immune diseases reveals distinct gene associations that converge on common pathways.

Genome-wide association studies (GWAS) have mapped thousands of susceptibility loci associated with immune-mediated diseases. To assess the extent of the genetic sharing across nine immune-mediated diseases we apply genomic structural equation modelling to GWAS data from European populations. We identify three disease groups: gastrointestinal tract diseases, rheumatic and systemic diseases, and allergic diseases. Although loci associated with the disease groups are highly specific, they converge on perturbing the same pathways. Finally, we test for colocalization between loci and single-cell eQTLs derived from peripheral blood mononuclear cells. We identify the causal route by which 46 loci predispose to three disease groups and find evidence for eight genes being candidates for drug repurposing. Taken together, here we show that different constellations of diseases have distinct patterns of genetic associations, but that associated loci converge on perturbing different nodes in T cell activation and signalling pathways.

Participation bias in the UK Biobank distorts genetic associations and downstream analyses.

While volunteer-based studies such as the UK Biobank have become the cornerstone of genetic epidemiology, the participating individuals are rarely representative of their target population. To evaluate the impact of selective participation, here we derived UK Biobank participation probabilities on the basis of 14 variables harmonized across the UK Biobank and a representative sample. We then conducted weighted genome-wide association analyses on 19 traits. Comparing the output from weighted genome-wide association analyses (neffective = 94,643 to 102,215) with that from standard genome-wide association analyses (n = 263,464 to 283,749), we found that increasing representativeness led to changes in SNP effect sizes and identified novel SNP associations for 12 traits. While heritability estimates were less impacted by weighting (maximum change in h2, 5%), we found substantial discrepancies for genetic correlations (maximum change in rg, 0.31) and Mendelian randomization estimates (maximum change in ßSTD, 0.15) for socio-behavioural traits. We urge the field to increase representativeness in biobank samples, especially when studying genetic correlates of behaviour, lifestyles and social outcomes.

An atlas of genetic scores to predict multi-omic traits.

The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables highly cost-effective and powerful analyses for studies that do not have multi-omics1. Here we examine a large cohort (the INTERVAL study2; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank3 to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.

Mendelian randomisation identifies priority groups for prophylactic EBV vaccination.

BACKGROUND: Epstein Barr virus (EBV) infects ~ 95% of the population worldwide and is known to cause adverse health outcomes such as Hodgkin's, non-Hodgkin's lymphomas, and multiple sclerosis. There is substantial interest and investment in developing infection-preventing vaccines for EBV. To effectively deploy such vaccines, it is vital that we understand the risk factors for infection. Why particular individuals do not become infected is currently unknown. The current literature, describes complex, often conflicting webs of intersecting factors-sociodemographic, clinical, genetic, environmental-, rendering causality difficult to decipher. We aimed to use Mendelian randomization (MR) to overcome the issues posed by confounding and reverse causality to determine the causal risk factors for the acquisition of EBV. METHODS: We mapped the complex evidence from the literature prior to this study factors associated with EBV serostatus (as a proxy for infection) into a causal diagram to determine putative risk factors for our study. Using data from the UK Biobank of 8422 individuals genomically deemed to be of white British ancestry between the ages of 40 and 69 at recruitment between the years 2006 and 2010, we performed a genome wide association study (GWAS) of EBV serostatus, followed by a Two Sample MR to determine which putative risk factors were causal. RESULTS: Our GWAS identified two novel loci associated with EBV serostatus. In MR analyses, we confirmed shorter time in education, an increase in number of sexual partners, and a lower age of smoking commencement, to be causal risk factors for EBV serostatus. CONCLUSIONS: Given the current interest and likelihood of a future EBV vaccine, these factors can inform vaccine development and deployment strategies by completing the puzzle of causality. Knowing these risk factors allows identification of those most likely to acquire EBV, giving insight into what age to vaccinate and who to prioritise when a vaccine is introduced.

Impact of cultural and genetic structure on food choices along the Silk Road.

The complex interplay between genetics, culture, and environment forms an individual's biology, influencing their behavior, choices, and health. However, to what extent information derived from this intertwined network could be quantitatively summarized to provide a glance at an individual's lifestyle is difficult to say. Here, we focused on dietary preferences as cultural proxies and genome-wide data of 543 individuals from six historical Silk Road countries: Georgia, Armenia, Azerbaijan, Uzbekistan, Kazakhstan, and Tajikistan. These lands favored the dispersal of innovations, foods, and DNA halfway across Eurasia, thus representing an ideal subject to explore interactions of cultural factors and genetic ancestry. We used discriminant analysis of principal components to infer cultural clusters, where mixed memberships are allowed. Five different clusters emerged. Of these, clusters 1 and 3, driven by aversion to pork and alcoholic beverages, mirrored genetic admixture patterns with the exception of Azerbaijan, which shares preferences supported by Islamic culture with Eastern countries. Cluster 3 was driven by protein-rich foods, whose preference was significantly related to steppe pastoralist ancestry. Sex and age were secondary clustering factors, with clusters formed by male and young individuals being related to alcohol preference and a reduced liking for vegetables. The soft clustering approach enabled us to model and summarize the individual's dietary information in short and informative vectors, which show meaningful interaction with other nondietary attributes of the studied individuals. Encoding other cultural variables would help summarize an individual's culture quantitatively, thus ultimately supporting its inclusion as a covariate in future association studies.

Using genetic variation to disentangle the complex relationship between food intake and health outcomes.

Diet is considered as one of the most important modifiable factors influencing human health, but efforts to identify foods or dietary patterns associated with health outcomes often suffer from biases, confounding, and reverse causation. Applying Mendelian randomization in this context may provide evidence to strengthen causality in nutrition research. To this end, we first identified 283 genetic markers associated with dietary intake in 445,779 UK Biobank participants. We then converted these associations into direct genetic effects on food exposures by adjusting them for effects mediated via other traits. The SNPs which did not show evidence of mediation were then used for MR, assessing the association between genetically predicted food choices and other risk factors, health outcomes. We show that using all associated SNPs without omitting those which show evidence of mediation, leads to biases in downstream analyses (genetic correlations, causal inference), similar to those present in observational studies. However, MR analyses using SNPs which have only a direct effect on the exposure on food exposures provided unequivocal evidence of causal associations between specific eating patterns and obesity, blood lipid status, and several other risk factors and health outcomes.

Large-scale GWAS of food liking reveals genetic determinants and genetic correlations with distinct neurophysiological traits.

We present the results of a GWAS of food liking conducted on 161,625 participants from the UK-Biobank. Liking was assessed over 139 specific foods using a 9-point scale. Genetic correlations coupled with structural equation modelling identified a multi-level hierarchical map of food-liking with three main dimensions: "Highly-palatable", "Acquired" and "Low-caloric". The Highly-palatable dimension is genetically uncorrelated from the other two, suggesting that independent processes underlie liking high reward foods. This is confirmed by genetic correlations with MRI brain traits which show with distinct associations. Comparison with the corresponding food consumption traits shows a high genetic correlation, while liking exhibits twice the heritability. GWAS analysis identified 1,401 significant food-liking associations which showed substantial agreement in the direction of effects with 11 independent cohorts. In conclusion, we created a comprehensive map of the genetic determinants and associated neurophysiological factors of food-liking.

Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Genome-wide Association Study of Liking for Several Types of Physical Activity in the UK Biobank and Two Replication Cohorts.

INTRODUCTION: A lack of physical activity (PA) is one of the most pressing health issues today. Our individual propensity for PA is influenced by genetic factors. Stated liking of different PA types may help capture additional and informative dimensions of PA behavior genetics. METHODS: In over 157,000 individuals from the UK Biobank, we performed genome-wide association studies of five items assessing the liking of different PA types, plus an additional derived trait of overall PA-liking. We attempted to replicate significant associations in the Netherlands Twin Register (NTR) and TwinsUK. Additionally, polygenic scores (PGS) were trained in the UK Biobank for each PA-liking item and for self-reported PA behavior, and tested for association with PA in the NTR. RESULTS: We identified a total of 19 unique significant loci across all five PA-liking items and the overall PA-liking trait, and these showed strong directional consistency in the replication cohorts. Four of these loci were previously identified for PA behavior, including CADM2 , which was associated with three PA-liking items. The PA-liking items were genetically correlated with self-reported ( rg = 0.38-0.80) and accelerometer ( rg = 0.26-0.49) PA measures, and with a wide range of health-related traits. Each PA-liking PGS significantly predicted the same PA-liking item in NTR. The PGS of liking for going to the gym predicted PA behavior in the NTR ( r2 = 0.40%) nearly as well as a PGS based on self-reported PA behavior ( r2 = 0.42%). Combining the two PGS into a single model increased the r2 to 0.59%, suggesting that PA-liking captures distinct and relevant dimensions of PA behavior. CONCLUSIONS: We have identified the first loci associated with PA-liking and extended our understanding of the genetic basis of PA behavior.

Identification of genetic loci simultaneously associated with multiple cardiometabolic traits.

BACKGROUND AND AIMS: Cardiometabolic disorders (CMD) arise from a constellation of features such as increased adiposity, hyperlipidemia, hypertension and compromised glucose control. Many genetic loci have shown associations with individual CMD-related traits, but no investigations have focused on simultaneously identifying loci showing associations across all domains. We therefore sought to identify loci associated with risk across seven continuous CMD-related traits. METHODS AND RESULTS: We conducted separate genome-wide association studies (GWAS) for systolic and diastolic blood pressure (SBP/DBP), hemoglobin A1c (HbA1c), low- and high- density lipoprotein cholesterol (LDL-C/HDL-C), waist-to-hip-ratio (WHR), and triglycerides (TGs) in the UK Biobank (N = 356,574-456,823). Multiple loci reached genome-wide levels of significance (N = 145-333) for each trait, but only four loci (in/near VEGFA, GRB14-COBLL1, KLF14, and RGS19-OPRL1) were associated with risk across all seven traits (P < 5 × 10-8). We sought replication of these four loci in an independent set of seven trait-specific GWAS meta-analyses. GRB14-COBLL1 showed the most consistent replication, revealing nominally significant associations (P < 0.05) with all traits except DBP. CONCLUSIONS: Our analyses suggest that very few loci are associated in the same direction of risk with traits representing the full spectrum of CMD features. We identified four such loci, and an understanding of the pathways between these loci and CMD risk may eventually identify factors that can be used to identify pathologic disturbances that represent broadly beneficial therapeutic targets.

Mendelian Randomization Identifies the Potential Causal Impact of Dietary Patterns on Circulating Blood Metabolites.

Nutrition plays an important role in the development and progress of several health conditions, but the exact mechanism is often still unclear. Blood metabolites are likely candidates to be mediating these relationships, as their levels are strongly dependent on the frequency of consumption of several foods/drinks. Understanding the causal effect of food on metabolites is thus of extreme importance. To establish these effects, we utilized two-sample Mendelian randomization using the genetic variants associated with dietary traits as instrumental variables. The estimates of single-nucleotide polymorphisms' effects on exposures were obtained from a recent genome-wide association study (GWAS) of 25 individual and 15 principal-component dietary traits, whereas the ones for outcomes were obtained from a GWAS of 123 blood metabolites measured by nuclear magnetic resonance spectroscopy. We identified 413 potentially causal links between food and metabolites, replicating previous findings, such as the association between increased oily fish consumption and higher DHA, and highlighting several novel associations. Most of the associations were related to very-low-density, intermediate-density (IDL), and low-density lipoproteins (LDL). For example, we found that constituents of IDL particles and large LDL particles were raised by coffee and alcohol while lowered by an overall healthier diet and fruit consumption. Our findings provide a strong base of evidence for planning future RCTs aimed at understanding the role of diet in determining blood metabolite levels.

Genetic analyses identify widespread sex-differential participation bias.

Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of participation bias is challenging since it requires the genotypes of unseen individuals. Here we demonstrate that it is possible to estimate comparative biases by performing a genome-wide association study contrasting one subgroup versus another. For example, we showed that sex exhibits artifactual autosomal heritability in the presence of sex-differential participation bias. By performing a genome-wide association study of sex in approximately 3.3 million males and females, we identified over 158 autosomal loci spuriously associated with sex and highlighted complex traits underpinning differences in study participation between the sexes. For example, the body mass index-increasing allele at FTO was observed at higher frequency in males compared to females (odds ratio = 1.02, P = 4.4 × 10-36). Finally, we demonstrated how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

Physician-Confirmed and Administrative Definitions of Stroke in UK Biobank Reflect the Same Underlying Genetic Trait.

BACKGROUND: Stroke in UK Biobank (UKB) is ascertained via linkages to coded administrative datasets and self-report. We studied the accuracy of these codes using genetic validation. METHODS: We compiled stroke-specific and broad cerebrovascular disease (CVD) code lists (Read V2/V3, ICD-9/-10) for medical settings (hospital, death record, primary care) and self-report. Among 408,210 UKB participants, we identified all with a relevant code, creating 12 stroke definitions based on the code type and source. We performed genome-wide association studies (GWASs) for each definition, comparing summary results against the largest published stroke GWAS (MEGASTROKE), assessing genetic correlations, and replicating 32 stroke-associated loci. RESULTS: The stroke case numbers identified varied widely from 3,976 (primary care stroke-specific codes) to 19,449 (all codes, all sources). All 12 UKB stroke definitions were significantly correlated with the MEGASTROKE summary GWAS results (rg.81-1) and each other (rg.4-1). However, Bonferroni-corrected confidence intervals were wide, suggesting limited precision of some results. Six previously reported stroke-associated loci were replicated using ≥1 UKB stroke definition. CONCLUSIONS: Stroke case numbers in UKB depend on the code source and type used, with a 5-fold difference in the maximum case-sample size. All stroke definitions are significantly genetically correlated with the largest stroke GWAS to date.

More Than Smell-COVID-19 Is Associated With Severe Impairment of Smell, Taste, and Chemesthesis.

Recent anecdotal and scientific reports have provided evidence of a link between COVID-19 and chemosensory impairments, such as anosmia. However, these reports have downplayed or failed to distinguish potential effects on taste, ignored chemesthesis, and generally lacked quantitative measurements. Here, we report the development, implementation, and initial results of a multilingual, international questionnaire to assess self-reported quantity and quality of perception in 3 distinct chemosensory modalities (smell, taste, and chemesthesis) before and during COVID-19. In the first 11 days after questionnaire launch, 4039 participants (2913 women, 1118 men, and 8 others, aged 19-79) reported a COVID-19 diagnosis either via laboratory tests or clinical assessment. Importantly, smell, taste, and chemesthetic function were each significantly reduced compared to their status before the disease. Difference scores (maximum possible change ±100) revealed a mean reduction of smell (-79.7 ± 28.7, mean ± standard deviation), taste (-69.0 ± 32.6), and chemesthetic (-37.3 ± 36.2) function during COVID-19. Qualitative changes in olfactory ability (parosmia and phantosmia) were relatively rare and correlated with smell loss. Importantly, perceived nasal obstruction did not account for smell loss. Furthermore, chemosensory impairments were similar between participants in the laboratory test and clinical assessment groups. These results show that COVID-19-associated chemosensory impairment is not limited to smell but also affects taste and chemesthesis. The multimodal impact of COVID-19 and the lack of perceived nasal obstruction suggest that severe acute respiratory syndrome coronavirus strain 2 (SARS-CoV-2) infection may disrupt sensory-neural mechanisms.

Variants associated with HHIP expression have sex-differential effects on lung function.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 ß=0.028 [SE 0.0022] litres) than females (ß=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Coffee Consumption and Kidney Function: A Mendelian Randomization Study.

RATIONALE & OBJECTIVE: Chronic kidney disease (CKD) is a leading cause of morbidity and mortality worldwide, with limited strategies for prevention and treatment. Coffee is a complex mixture of chemicals, and consumption has been associated with mostly beneficial health outcomes. This work aimed to determine the impact of coffee consumption on kidney function. STUDY DESIGN: Genome-wide association study (GWAS) and Mendelian randomization. SETTING & PARTICIPANTS: UK Biobank baseline data were used for a coffee consumption GWAS and included 227,666 participants. CKDGen Consortium data were used for kidney outcomes and included 133,814 participants (12,385 cases of CKD) of mostly European ancestry across various countries. EXPOSURE: Coffee consumption. OUTCOMES: Estimated glomerular filtration rate (eGFR), CKD GFR categories 3 to 5 (G3-G5; eGFR<60mL/min/1.73m2), and albuminuria. ANALYTICAL APPROACH: GWAS to identify single-nucleotide polymorphisms (SNPs) associated with coffee consumption in UK Biobank and use of those SNPs in Mendelian randomization analyses of coffee consumption and kidney outcomes in CKDGen. RESULTS: 2,126 SNPs were associated with coffee consumption (P<5×10-8), 25 of which were independent and available in CKDGen. Drinking an extra cup of coffee per day conferred a protective effect against CKD G3-G5 (OR, 0.84; 95% CI, 0.72-0.98; P=0.03) and albuminuria (OR, 0.81; 95% CI, 0.67-0.97; P=0.02). An extra cup was also associated with higher eGFR (ß=0.022; P=1.6×10-6) after removal of 3 SNPs responsible for significant heterogeneity (Cochran Q P = 3.5×10-15). LIMITATIONS: Assays used to measure creatinine and albumin varied between studies that contributed data and a sex-specific definition was used for albuminuria rather than KDIGO guideline recommendations. CONCLUSIONS: This study provides evidence of a beneficial effect of coffee on kidney function. Given widespread coffee consumption and limited interventions to prevent CKD incidence and progression, this could have significant implications for global public health in view of the increasing burden of CKD worldwide.

Associations of autozygosity with a broad range of human phenotypes.

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

A genome-wide association study identifies an association between variants in EFCAB4B gene and periodontal disease in an Italian isolated population.

BACKGROUND AND OBJECTIVE: Periodontitis in one of the most prevalent dental diseases. Despite numerous studies have investigated its aetiopathogenetic factors, few works have focused on its genetic predisposition and most of them took into account only candidate genes. Therefore, we conducted a Genome Wide Association Study in an Italian isolated population aimed at uncovering genetic variants that predispose to this disorder. METHODS: Diagnosis of chronic periodontitis was made following the criteria of the American Academy of Periodontology. Patients with chronic periodontitis were grouped into different categories: slight, severe, localized and generalized. A control group composed by people without signs of periodontitis or gingivitis was defined. DNA was genotyped using 370k Illumina chips. Linear mixed model regression was used to test the association between each single nucleotide polymorphism (SNP) (independent variable) and the periodontitis status (dependent variable), controlling for confounders sex, age and smoking. The genomic kinship matrix was also used as random effect. RESULTS: Four SNPs on the gene EFCAB4B resulted significantly associated to localized periodontitis (P < 5 × 10-8 ), with the best hit on the rs242016 SNP (P = 1.5 × 10-8 ). CONCLUSION: We have identified a novel significant association between the EFCAB4B gene and localized periodontitis. These results open a new perspective in the understanding of genetic factors contributing to this common disorder.